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Abstract 

We present studies of two-body and three-body charmed baryonic B decays in a sample of 232 
million BB pairs collected with the BABAR detector at the PEP-II e + e~ storage ring. The branching 
fractions of the decays B° — > A^p and B~ — > A^fm~ are measured to be (2.15±0.36±0.13±0.56) x 
10" 5 and (3.53 ± 0.18 ± 0.31 ± 0.92) x 10" 4 , respectively. The uncertainties quoted are statistical, 
systematic, and from the — * pK~ir + branching fraction. We observe a baryon-antibaryon 
threshold enhancement in the A+p invariant mass spectrum of the three-body mode and measure 
the ratio of the branching fractions to be B(B~ — > A+pir~)/B(B° — >• A+p) = 16.4±2.9± 1.4. These 
results are preliminary. 
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1 INTRODUCTION 



Charmed baryonic B decays are experimentally accessible and provide a way to check predictions 
given by various theoretical models for exclusive baryonic B decays. There is theoretical inter- 
est in the suppression of the two-body baryonic decay rates compared to three-body decay rates 
and the possible connection to production mechanisms for baryons in B decays. Analysis of the 
charmed three-body baryonic B decay reveals that the invariant mass of the baryon-antibaryon 
system is peaked near threshold [1]. Charmless two-body baryonic B decays (which have not 
yet been observed [2,3]) may be used to measure direct CP violation in the B system. Their 
charmed counterparts, however, have branching fractions at least an order of magnitude higher 
than the charmless modes, and thus can help distinguish between theoretical models that predict 
the charmless decay rates of B mesons to baryons. The Feynman diagrams for these decays are 
shown in Figure 1, in which the B meson decays weakly via internal W emission to A~£p(ir). 

Charmed baryonic B decays have recently been measured by the CLEO [4] and Belle [1,5,6] 
Collaborations. In particular, the Belle Collaboration has measured the branching fractions of the 
modes 5 B° — > A+p (using 85 million BB pairs) [5] and B~ — > A~£p7t~ (using 152 million BB 
pairs) [1]: 

B(B° -► A+p) = (2.19±g;tg ± 0.32 ± 0.57) x 1(T 5 and 

B(B- -► A+pn-) = (20.1 ± 1.5 ± 2.0 ± 5.2) x 10" 5 , 

where the errors are statistical, systematic, and from the — ► pK~ir + branching fraction, respec- 
tively. BABAR has collected nearly three times the data used in the Belle analysis of the two-body 
mode, and we can therefore perform a more precise measurement of this branching fraction. For 
now, the measurement errors are dominated by the 26% fractional error on B(A+ — > pK~ir + ) = 
(5.0 ± 1.3)% [7], but this uncertainty cancels in the ratio of the three-body to two-body branching 
fractions. 

The excess of events near the baryon-antibaryon production threshold seen by Belle in B~ — > 
A+pir~ has also been observed in B° — ► Apir~ [8] and several B — > ppX [9, 10] modes. In refer- 

5 Throughout this paper, whenever a mode is given, the charge conjugate is also implied. 
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ence [11] a qualitative explanation of the larger three-body branching fraction in conjunction with 
this threshold effect is given. In the two-body decay, the invariant mass of the baryon-antibaryon 
is simply rag, whereas in the three-body decay, the invariant mass of the baryon-antibaryon can 
be lower, allowing the baryon-antibaryon to form a quasi-resonance near threshold. The third 
daughter, the meson, carries away much of the energy. The result, regardless of the interpretation 
of the threshold enhancement, is that the B favors three-body baryonic decay modes by an order 
of magnitude over two-body modes. 

In this analysis, we measure the branching fractions for B° — > A+p and B~ — ► Afpir~ and 
observe the threshold enhancement in the baryon-antibaryon system of the the B~ — > A+pir~ 
mode. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II e + e~ stor- 
age ring. The data sample used comprises an integrated luminosity of 210 fb -1 (232 million BB 
pairs) collected from e + e~ collisions at the T(4S) resonance. The BABAR detector is described else- 
where [12]. Exclusive B meson decays are simulated with the Monte Carlo (MC) event generator 
EvtGen [13] and hadronization (e.g. for continuum qq events) is simulated with Jetset7.4 [14]. 
The detector is modeled using the GEANT4 simulation package [15]. 

3 ANALYSIS METHOD 

3.1 Candidate Selection 

We reconstruct A£ candidates in the decay mode A£ — > pK~ir + , applying a geometric constraint 
on the pK~7r + vertex, which is required to have a x 2 probability greater then 0.1%. The pK~ir + 
invariant mass must be between 2.275 and 2.295 GeV/c 2 . The pK~ir + candidates are constrained to 
the mass of the A+ [7], which provides better resolution in the kinematic variable AE = E* B — y/s/2, 
where E* B is the B candidate energy in the e + e~ center-of-mass (CM) frame and \/s is the total 
CM energy. Af candidates are then combined in a geometric fit with a p (and tt) to form a B° 
(B~) candidate for the B° -► A+p (B~ -► A+ pn ) mode. The \ 2 probability for the fit to the full 
decay tree must be greater than 0.1%. 

Daughter p, K, and ir candidates must be well-reconstructed in the drift chamber and are 
identified with likelihood-based particle selectors using information from the silicon vertex tracker, 
drift chamber, and ring-imaging Cerenkov detector. Several requirements differ between the two- 
and three-body modes. The pions in the B~ — > A+pir" mode have lower momenta; therefore, 
we apply looser drift chamber tracking requirements to improve the efficiency in several areas of 
the A+pw~ Dalitz plane. The daughter particles in both decay modes have very loose particle 
identification requirements with two exceptions: 1) the pions in B~ — ► A+pm~ are required to 
satisfy stronger kaon and electron rejection criteria, and 2) the B daughter p in B° — ► A+p must 
pass a tight constraint on the likelihood that the track is a proton and stronger electron rejection. 

We construct a linear (Fisher) discriminant T from several event-shape variables to provide 
continuum suppression: | cos 9*\ (9* is the angle of the B candidate momentum vector with respect 
to the beam axis in the e + e~ CM frame), | cos#^ r | (9f hr is the angle of the B candidate thrust 
axis with respect to the beam axis in the e + e~ CM frame), and the summed momentum of the 
rest of the charged and neutral particles in the event in nine cones centered around the thrust 
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axis of the B candidate. The requirement on T provides powerful background rejection (72.8%) 
for the two-body mode, but is less effective for the three-body mode (28.0%) due to a larger 
component of combinatoric B backgrounds compared to the continuum component. These values 
were determined by maximizing N s /y/N s + iV&, where N s is the number of signal events based on 
signal MC samples and TV;, is the number of background events in AE upper sidebands (AE > 
0.1 GeV and 5.2 < m E s < 5.29GeV/c 2 ) in data. 

We identify signal candidates using AE and the beam-energy-substituted mass 

mEs = \J ((s/2 + Pi • PB) 2 /Ef — p 2 B ), where (-E^pj) is the four-momentum of the initial e + e~ 
system and ps is the momentum of the B candidate, both measured in the laboratory frame. The 
distribution of AE vs. rriEs f° r both modes is shown in Figure 2. We define the fit region to be 
-0.1 < AE < 0.1 GeV and 5.20 < m E s < 5.29 GeV/ c 2 (also indicated in Figure 2). This excludes 
the AE sideband used in the optimization and the region below —0.1 GeV in AE, which contains 
backgrounds that peak in mES but are shifted in AE. These backgrounds are from B — > A c pTrir 
(B — ► A c p7r) events where a B daughter n is not included in the B candidate, mimicking the mode 
of interest: B~ — ► A+pir~ (B° — > A+p). Studies of exclusive MC samples of these backgrounds 
indicate that much of the contribution is from B — > E c pit (B — > E c p) where a 7r° or slow charged ir 
from the S c — > A c ir decay is missed. MC samples comprised of continuum qq events and B meson 
decays were studied to rule out any background that peaks in both AE and mES- 
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(a) B° candidates in data for B° — ► yl+p 



(b) B candidates in data for B — » A+pir 



Figure 2: Distribution of AE vs. mss of candidates in data for both the two-body (a) and 
three-body (b) decay modes. The fit regions are indicated. 



3.2 B° — ► A+p Maximum Likelihood Fit 

In the analysis of the two-body B° — > yl+p mode, we did not look at the signal region until the 
event selection criteria and fit procedures were determined. The efficiency for reconstructing and 
selecting B° — > A^p candidates is 20.2%, and is determined from a fit to the B° — ► yl^p signal 
MC sample. A 2-D unbinned maximum likelihood fit is performed in AE and mES to extract the 
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Figure 3: Projections of the 2-D fit in AE and toes f° r ^tP candidates satisfying \AE\ < 0.04 GeV 
(left) and toes > 5.27 GeV/ c 2 (right). The signal yield is 50 ± 8 events, with a significance of 9Aa. 



number of signal events. The background is described by the product of a linear function in AE 
and a threshold function [16] in toes! the signal is described by a single Gaussian distribution in 
each dimension. All parameters except the toes threshold are unconstrained in the fit to data. 
We validate the fitting procedure on a combined sample of signal MC events (over a range of the 
expected number of signal events) and "toy" MC events (generated according to the shape of the 
continuum and BB MC background events) to ensure that the fit is robust and unbiased. 

The results of the fit to data are shown in Figure 3; we obtain 50 ± 8 signal events and a 
significance of s/—2 In (Co/C max ) = 9.4<r, where C max is the maximum likelihood from the fit 
result and Cq is the maximum likelihood when the signal yield is fixed to zero. The mean in AE 
is shifted slightly below zero (—4.2 ± 2.7 MeV); this shift is in the appropriate direction given that 
the mass is constrained to the 2004 PDG value [7] which is approximately 1.5 MeV lower than 
the most recent measurement [17]. The AE resolution, 15.4 ± 2.1 MeV, is slightly larger than, but 
consistent with, the resolution in MC (13.6 ± 0.1 MeV). 

3.3 B~ — > A+ptt~ Maximum Likelihood Fit 

For the three-body B~ — > A^pn~ mode, a 2-D unbinned maximum likelihood fit is also performed. 
Again, all parameters except the toes threshold are unconstrained in the fit to data. The back- 
ground PDF is the same as in the two-body mode, but the signal PDF consists of a Gaussian in 
AE times a Gaussian in toeSj where a correlation is allowed between the two observables. This 
was not necessary in the two-body mode due to the limited number of signal events. The signal 
PDF also contains an additional uncorrelated Gaussian component in AE with the same mean as 
the correlated Gaussian but an independent width. This signal PDF was chosen from a study of 
B~ — > A+pir~ signal MC events along with extensive studies of various PDFs using a combined 
sample of signal MC and toy MC events. These studies showed this PDF to have the smallest bias: 
—8 ± 2 events for 500 total signal events (the level of bias is consistent for a range of signal events). 
The result of the fit to data with this PDF is shown in Figure 4. The signal yield from the fit is 
571 ± 34 events and the AE resolution (RMS) is 19 ± 3 MeV. 
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Figure 4: Projections of the 2-D fit in tttes and AE, for A+pir candidates satisfying \AE\ < 
0.030 GeV (left) and m ES > 5.27GeV/c 2 (right). This 2-D fit is used to extract the likelihood that 
each event is a signal or background event. The signal yield is 571 ± 34 events. 



3.4 B — > A+pir Yield Extraction and Efficiency Correction 

We use the a Vlot technique [18] (a sophisticated background subtraction method) to project out the 
signal and background distributions separately based on the 2-D fit to AE and tttes- We calculate 
a signal weight for each event i according to the following equation: 

w _ fs(m E Si, AEj) + Vsbfb(m E Si, AEj) 
1 ~ NsfsirriEs^AEi) + N b f b (m ESi , AE*) ' 1 ' 

where W% is the s Vlot weight, N s (N b ) is the number of fitted signal (background) events, and 
fs(fb) is the signal (background) PDF. Y sb is the off-diagonal element of a 2 x 2 covariance matrix 
calculated directly from data, with all parameters fixed to their fitted values except for the signal 
and background yields. A background weight for each event can be calculated in an analogous 
manner. The result of this method is that each event is assigned a signal and background weight, 
which can be plotted for any quantity that is uncorrelated with AE and ttt.es- The quantities 
of interest that satisfy this requirement are the invariant masses m^d ■> where d{ is any of the B 
daughters A£, p, tt~ . The correlations of AE and tttes with these quantities are less than 5%. The 
s Vlot method relies on using the events in the entire fit region to provide good sampling of both 
signal and background. However, (background) events that have an invariant A^pir^ mass far from 
the mass of the B meson have a different kinematically allowed Dalitz region than (signal) events 
with an invariant A^pn~ mass close to m/3. We calculate m^d with a B mass constraint so that 
all of the B candidates in the fit region lie in the same Dalitz region. 

The detection efficiency for B~ — ► A^pir~ events varies significantly across the Dalitz plane. 
Therefore, using the average nonresonant MC efficiency (15.3%) to calculate the branching fraction 
for this mode is insufficient. Instead, an efficiency correction is applied to each signal event based 
on its location in the Dalitz plane. We divide the physical region into 215 equal-size bins and 
determine the efficiency in each bin; a plot of this efficiency for vs. nt^ c7r is shown in Figure 5. 
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< c n (GeV 2 /c 4 ) 

Figure 5: Binned efficiency for vs. rn\ cn in the kinematically allowed region of the Dalitz 
plane. 



There are noticeable deficiencies in the lower left (right) corners of the A+p~Tr~ Dalitz plane, where 
the 7r (A c ) candidates have low momentum in the B rest frame. The looser tracking requirements 
on the pions help to compensate for this effect, but do not eliminate it entirely. We build on the 
/Plot formalism to individually correct each event by an additional weight, l/e a p, where e a p is 
the efficiency in bin (a, (3). We define an "effective" efficiency (e e ff) as the signal yield from the fit 
divided by the number of s'P/ot-weighted, efficiency-corrected events. The effective efficiency for 
selecting B candidates in the three-body mode is e e s = 14.2%. 

4 SYSTEMATIC STUDIES 

Various sources of systematic uncertainties have been investigated, including those related to the 
total number of BB pairs in data, the method used to determine the efficiency from MC, and the 
fitting procedures. These are summarized for both modes in Table 1. Note that for the branching 
fraction measurement of B° — ► Afp, the statistical error dominates over the total systematic 
uncertainty. 

The systematic uncertainty on the number of BB pairs produced by BABAR is 1.1%. 

There are several sources of systematic uncertainty related to the efficiency determinations. 
The statistical uncertainty due to the number of signal MC events contributes a 1.0% systematic 
error on the efficiency. Tracking efficiency systematic errors are based on studies of r decays, 
which yield an uncertainty of 0.8% per track. However, this is reduced to 0.6% for the higher 
momentum B daughter p in the two-body mode and increased to 1.4% for the lower momentum 
pions in the three-body mode. Particle identification is determined using large control samples, 
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Table 1: Summary of the contributions to the total systematic uncertainty. The total is determined 
by adding the uncertainty from each source in quadrature. 



Source 


Systematic Uncertainty 




B° -+ A+p 


B -»■ A+pn~ 




1.1% 


1.1% 


MC statistics 


1.0% 


4.7% 


Dalitz binning 




2.0% 


Tracking 


3.0% 


5.2% 


PID 


4.7% 


1.2% 


Fitting 


2.2% 


4.5% 


Total 


6.0% 


8.7% 



which may differ from the modes we are investigating due to the higher multiplicities of these 
charmed baryonic B decays and other subtleties. Differences between the momentum spectra and 
angular distributions of the daughter particles compared to those in the control samples are used 
to assess a systematic uncertainty on the efficiency due to particle identification. In B° — ► A+p, we 
assign a 2.5% systematic uncertainty to the B daughter p, and a (1.5 to 1.7)% uncertainty on the 
other daughter particles. In B~ — > A+pir~ , the systematic uncertainty due to particle identification 
varies from (0.1 to 0.9)% per track; the total is 1.2%. 

The fitting systematics are studied by varying the background shape in AE and the endpoint 
of the threshold function in ttjes- This yields a systematic uncertainty of 2.2% for the two-body 
mode and 0.9% for the three-body mode. The branching fraction measurement of B~ — > A+pir~ 
has additional systematic uncertainties due to the signal PDF. We assign 4.3% due to the fit bias, 
the source of which is mostly in the the tails of AE. This systematic uncertainty compensates for 
the inability of the MC to accurately simulate the behavior of the events in these tails. For the 
three-body mode, we perform the fit to data with and without a correlation between AE and tbes 
in the signal PDF, yielding a systematic uncertainty of 1.1%. 

5 RESULTS 

The branching fraction of B° —> A+p [19], measured with the sample of 232 million BB pairs, is 

B(B° A+p) = (2.15 ± 0.36 ± 0.13 ± 0.56) x 10" 5 , 

where the errors are statistical, systematic, and from the A+ — ► pK~ir+ branching fraction, respec- 
tively. The significance of the signal is 9.4<r. This measurement is consistent with a previous Belle 
measurement of (2.19:^49 ± °- 32 ± °- 57 ) x 10 ~ 5 made with 85 mil li° n BB pairs. The systematic 
uncertainty is much lower (6% compared to 15%) than that for the Belle measurement. We also 
find this measurement to be consistent with the predicted limit from reference [20]: B(B° — > A+p) 
< 7.9 x 10" 6 \g/5\ 2 , where \g\ = 6-10. 

We calculate the total branching fraction of B~ — > A+pit~ as follows: 
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Figure 6: mZ^ vs. rn 2 AclT for signal B candidates with s Tlot and efficiency correction weights applied. 
Bins with negative population are suppressed. 



B(B- A c V-,« = ——^-^ (2) 
= (3.53 ± 0.18 ± 0.31 ± 0.92) x 10~ 4 , 

where the fit bias b is 1.6%, Wi is the signal s Vlot weight and is the efficiency for event i, and N BB 
is the number of BB pairs. The uncertainties are statistical, systematic, and the error on the A+ -> 
pK~7r + branching fraction, respectively. This measurement is 3.5 a higher (assuming Gaussian 
statistics) than the Belle measurement of B(B~ -> A+p-K-) to t = (2.01 ± 0.15 ± 0.20 ± 0.52) x 10~ 4 . 
An examination of the Dalitz plot shows a systematic trend in that we measure consistently larger 
branching fractions in all regions. 

The Afpir~ Dalitz plane in data is shown in Figure 6 with s Vlot weights and efficiency cor- 
rections applied to each B candidate. We project this onto the va/y cV axis with the requirement 
mA c n > 2.6GeV/c 2 (to remove the contribution from the I7 C (2455)°) in Figure 7. We observe 
a baryon-antibaryon threshold enhancement in the A^p invariant mass spectrum, confirming the 
large body of evidence supporting the existence of these threshold enhancements in three-body 
baryonic B decays. 

We also report the ratio of the branching fractions of B~ — ► Afpn^ to B° — > A^p: 

^lT^%Ll = 16 . 4 ±2.9 ±1.4. (3) 
B(B° -» A+p) y ' 
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Figure 7: Dalitz plot projection onto the mj\ cP axis with the requirement my\ c7T > 2.6GeV/c 2 , 
removing the contribution from the I7 C (2455) . s Vlot weighted, efficiency-corrected signal events 
are shown. The baryon-antibaryon threshold enhancement is visible near 3.3GeV/c 2 . 



Table 2: Comparison of the yields, efficiencies (effective for the three-body decay), and branching 
fractions for B° — ► A+p and B~ — > Afp7r~. 



Mode 


Signal yield 


e (eS) 


B 


B° -> /L+p 


50 ±8 


20.2% 


(2.15 ± 0.36 ± 0.13 ± 0.56) x 10~ 5 


5~ — > A+pir~ 


571 ± 34 


14.2% 


(3.53 ± 0.18 ± 0.31 ± 0.92) x 10" 4 



The systematic uncertainties on the number of BB pairs, the Af daughter p and K tracking, and 
the A£ daughter K and B daughter p particle identification all cancel, as does the uncertainty on 
B(A^ — > pi^ _ 7r + ). This ratio is consistent with theoretical predictions. 

6 SUMMARY 

We report the branching fractions of two charmed baryonic B decay modes. Table 2 compares 
the yields, efficiencies, and branching fractions of the two modes. The total three-body branching 
fraction measured is significantly larger than that measured by Belle, but is still consistent with 
(and perhaps provides stronger evidence for) the observation that the three-body mode is enhanced 
over the two-body mode. The measurement of the ratio of three-body to two-body branching 
fractions and the observation of the baryon-antibaryon threshold enhancement aid in theoretical 
interpretations of baryon production in B decays. 
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